Dépliages et interprétation abstraite pour réseaux de régulation biologiques paramétrés

The analysis of dynamics of biological regulatory networks, notably signalling and gene regulatory networks, faces the uncertainty of the exact computational model. Indeed, most of the knowledge available concerns the existence of (possibly indirect) interactions between biological entities (species), e.g. proteins, RNAs, genes, etc. The details on how different regulators of a same target cooperate, and even more so on consistent rates for those interactions, however, are rarely available. In this regard, qualitative modelling approaches in the form of discrete regulatory networks, such as Boolean and Thomas networks, offer an appropriate level of abstraction for the biological regulatory network dynamics. As discrete regulatory networks are based on an influence graph, they require few additional parameters compared to classical quantitative models. Nevertheless, determining the discrete parameters is a well known challenge, and a major bottleneck for providing robust predictions from computational models.The influence graph of a regulatory network establishes dependencies for the evolution of each specie, specified by the directed edges of the graph. The dependencies alone, however, do not suffice to specify the logical function governing the evolution of a specie. Instead the logical functions associated to each specie, constrained by the influence graph, are encoded within the parameters of a discrete regulatory network. The space of admissible logical functions is then represented by a parametric regulatory network. On the one hand, parametric regulatory networks can be used for identification of parameter values for which the resulting discrete regulatory network satisfies given (dynamical) properties. Parameter identification of regulatory networks can thus be seen as a particular instance of model synthesis, in the constrained setting of the underlying influence graph. On the other hand, parametric regulatory networks may be analysed as a stand-alone model, for making predictions that are robust with respect to variability in the network.The analysis of parametric regulatory network dynamics is hampered by dual combinatorial explosion, of the state space and of the parameter space. In this thesis, we develop novel methods of parametric regulatory network analysis, in the form of specialised semantics, aimed at alleviating the combinatorial explosion. First, we introduce abstract interpretation for the set of admissible parameter evaluations (parametrisations).The abstraction allows us to represent any set of parametrisations by a constant size encoding, at the cost of a conservative over-approximation. Second, we lift partial order semantics in the form of unfolding from Petri nets to parametric regulatory networks. The influence graphs of biological regulatory networks tend to be relatively sparse, allowing for a lot of concurrency. This can be harnessed by partial order reduction methods to produce concise state space representations.The two approaches are aimed at tackling both aspects of the dual combinatorial explosion and are introduced in a compatible manner, allowing one to employ them simultaneously. Such application is supported by a prototype implementation used to conduct experiments on various parametric regulatory networks. We further consider refinements of the methods, such as an on-the-run model reduction method lifted to parametric regulatory networks from automata networks.L'analyse de la dynamique des réseaux de régulation biologique, notamment des réseaux de signalisation et de régulation génique, fait face à l'incertitude du modèle de calcul exact. En effet, la plupart des connaissances disponibles concernent l'existence d'interactions (éventuellement indirectes) entre des entités biologiques (espèces), par ex. protéines, ARN, gènes, etc. Les détails sur la manière dont les différents régulateurs d'une même cible coopèrent, et plus encore sur les taux cohérents pour ces interactions, sont cependant rarement disponibles. A cet égard, des approches de modélisation qualitative sous forme de réseaux de régulation discrets, tels que les réseaux booléens et Thomas, offrir un niveau d'abstraction approprié pour la dynamique du réseau de régulation biologique. Les réseaux de régulation discrets étant basés sur un graphe d'influence, ils nécessitent peu de paramètres supplémentaires par rapport aux modèles quantitatifs classiques. Néanmoins, la détermination des paramètres discrets est un défi bien connu et un goulot d'étranglement majeur pour fournir des prédictions robustes à partir de modèles informatiques.Le graphe d'influence d'un réseau de régulation établit des dépendances pour l'évolution de chaque espèce, spécifiées par les arêtes dirigées du graphe. Les dépendances seules, cependant, ne suffisent pas pour spécifier la fonction logique régissant l'évolution d'une espèce. Au lieu de cela, les fonctions logiques associées à chaque espèce, contraintes par le graphe d'influence, sont codées dans les paramètres d'un réseau de régulation discret. L'espace des fonctions logiques admissibles est alors représenté par un réseau de régulation paramétrique. D'une part, les réseaux de régulation paramétriques peuvent être utilisés pour l'identification de valeurs de paramètres pour lesquelles le réseau de régulation discret résultant satisfait des propriétés (dynamiques) données. L'identification des paramètres des réseaux de régulation peut ainsi être vue comme un exemple particulier de synthèse de modèle, dans le cadre contraint du graphe d'influence sous-jacent. D'autre part, les réseaux de régulation paramétriques peuvent être analysés comme un modèle autonome, pour faire des prédictions robustes vis-à-vis de la variabilité du réseau.L'analyse de la dynamique du réseau de régulation paramétrique est entravée par la double explosion combinatoire, de l'espace d'états et de l'espace des paramètres. Dans cette thèse, nous développons de nouvelles méthodes d'analyse de réseau de régulation paramétrique, sous forme de sémantique spécialisée, visant à atténuer l'explosion combinatoire. Tout d'abord, nous introduisons une interprétation abstraite de l'ensemble des évaluations de paramètres admissibles (paramétrisations). L'abstraction permet de représenter n'importe quel ensemble de paramétrisations par un encodage de taille constante, au prix d'une sur-approximation conservatrice. Deuxièmement, nous élevons la sémantique d'ordre partiel sous la forme d'un déploiement des réseaux de Petri vers des réseaux de régulation paramétriques. Les graphiques d'influence des réseaux de régulation biologique ont tendance à être relativement clairsemés, ce qui permet une grande concurrence. Cela peut être exploité par des méthodes de réduction d'ordre partiel pour produire des représentations d'espace d'état concises.Les deux approches visent à aborder les deux aspects de la double explosion combinatoire et sont introduites de manière compatible, ce qui permet de les utiliser simultanément. Une telle application est soutenue par une implémentation prototype utilisée pour mener des expériences sur divers réseaux de régulation paramétriques. Nous considérons en outre des raffinements des méthodes, comme une méthode de réduction de modèle à la volée portée aux réseaux de régulation paramétriques à partir de réseaux d'automates

Most Permissive Semantics of Boolean Networks

As shown in (http://dx.doi.org/10.1101/2020.03.22.998377), the usual update modes of Boolean networks (BNs), including synchronous and (generalized) asynchronous, fail to capture behaviors introduced by multivalued refinements. Thus, update modes do not allow a correct abstract reasoning on dynamics of biological systems, as they may lead to reject valid BN models.This technical report lists the main definitions and properties of the most permissive semantics of BNs introduced in http://dx.doi.org/10.1101/2020.03.22.998377. This semantics meets with a correct abstraction of any multivalued refinements, with any update mode. It subsumes all the usual updating modes, while enabling new behaviors achievable by more concrete models. Moreover, it appears that classical dynamical analyzes of reachability and attractors have a simpler computational complexity:- reachability can be assessed in a polynomial number of iterations. The computation of iterations is in NP in the very general case, and is linear when local functions are monotonic, or with some usual representations of functions of BNs (binary decision diagrams, Petri nets, automata networks, etc.). Thus, reachability is in P with locally-monotonic BNs, and P$^{\text{NP}}$ otherwise (instead of being PSPACE-complete with update modes);- deciding wherever a configuration belongs to an attractor is in coNP with locally-monotonic BNs, and coNP$^{\text{coNP}}$ otherwise (instead of PSPACE-complete with update modes).Furthermore, we demonstrate that the semantics completely captures any behavior achievable with any multilevel or ODE refinement of the BN; and the semantics is minimal with respect to this model refinement criteria: to any most permissive trajectory, there exists a multilevel refinement of the BN which can reproduce it.In brief, the most permissive semantics of BNs enables a correct abstract reasoning on dynamics of BNs, with a greater tractability than previously introduced update modes

Relational Differential Dynamic Logic

International audienceIn the field of quality assurance of hybrid systems (that combine continuous physical dynamics and discrete digital control), Platzer's differential dynamic logic (dL) is widely recognized as a deductive verification method with solid mathematical foundations and sophisticated tool support. Motivated by benchmarks provided by our industry partner , we study a relational extension of dL, aiming to formally prove statements such as "an earlier deployment of the emergency brake decreases the collision speed." A main technical challenge here is to relate two states of two dynamics at different time points. Our main contribution is a theory of suitable relational differential invariants (a relational extension of differential invariants that are central proof methods in dL), and a derived technique of time stretching. The latter features particularly high applicability, since the user does not have to synthesize a relational differential invariant out of the air. We derive new inference rules for dL from these notions, and demonstrate their use over a couple of automotive case studies

Dépliages et interprétation abstraite pour réseaux de régulation biologiques paramétrés

The analysis of dynamics of biological regulatory networks, notably signalling and gene regulatory networks, faces the uncertainty of the exact computational model. Indeed, most of the knowledge available concerns the existence of (possibly indirect) interactions between biological entities (species), e.g. proteins, RNAs, genes, etc. The details on how different regulators of a same target cooperate, and even more so on consistent rates for those interactions, however, are rarely available. In this regard, qualitative modelling approaches in the form of discrete regulatory networks, such as Boolean and Thomas networks, offer an appropriate level of abstraction for the biological regulatory network dynamics. As discrete regulatory networks are based on an influence graph, they require few additional parameters compared to classical quantitative models. Nevertheless, determining the discrete parameters is a well known challenge, and a major bottleneck for providing robust predictions from computational models.The influence graph of a regulatory network establishes dependencies for the evolution of each specie, specified by the directed edges of the graph. The dependencies alone, however, do not suffice to specify the logical function governing the evolution of a specie. Instead the logical functions associated to each specie, constrained by the influence graph, are encoded within the parameters of a discrete regulatory network. The space of admissible logical functions is then represented by a parametric regulatory network. On the one hand, parametric regulatory networks can be used for identification of parameter values for which the resulting discrete regulatory network satisfies given (dynamical) properties. Parameter identification of regulatory networks can thus be seen as a particular instance of model synthesis, in the constrained setting of the underlying influence graph. On the other hand, parametric regulatory networks may be analysed as a stand-alone model, for making predictions that are robust with respect to variability in the network.The analysis of parametric regulatory network dynamics is hampered by dual combinatorial explosion, of the state space and of the parameter space. In this thesis, we develop novel methods of parametric regulatory network analysis, in the form of specialised semantics, aimed at alleviating the combinatorial explosion. First, we introduce abstract interpretation for the set of admissible parameter evaluations (parametrisations).The abstraction allows us to represent any set of parametrisations by a constant size encoding, at the cost of a conservative over-approximation. Second, we lift partial order semantics in the form of unfolding from Petri nets to parametric regulatory networks. The influence graphs of biological regulatory networks tend to be relatively sparse, allowing for a lot of concurrency. This can be harnessed by partial order reduction methods to produce concise state space representations.The two approaches are aimed at tackling both aspects of the dual combinatorial explosion and are introduced in a compatible manner, allowing one to employ them simultaneously. Such application is supported by a prototype implementation used to conduct experiments on various parametric regulatory networks. We further consider refinements of the methods, such as an on-the-run model reduction method lifted to parametric regulatory networks from automata networks.L'analyse de la dynamique des réseaux de régulation biologique, notamment des réseaux de signalisation et de régulation génique, fait face à l'incertitude du modèle de calcul exact. En effet, la plupart des connaissances disponibles concernent l'existence d'interactions (éventuellement indirectes) entre des entités biologiques (espèces), par ex. protéines, ARN, gènes, etc. Les détails sur la manière dont les différents régulateurs d'une même cible coopèrent, et plus encore sur les taux cohérents pour ces interactions, sont cependant rarement disponibles. A cet égard, des approches de modélisation qualitative sous forme de réseaux de régulation discrets, tels que les réseaux booléens et Thomas, offrir un niveau d'abstraction approprié pour la dynamique du réseau de régulation biologique. Les réseaux de régulation discrets étant basés sur un graphe d'influence, ils nécessitent peu de paramètres supplémentaires par rapport aux modèles quantitatifs classiques. Néanmoins, la détermination des paramètres discrets est un défi bien connu et un goulot d'étranglement majeur pour fournir des prédictions robustes à partir de modèles informatiques.Le graphe d'influence d'un réseau de régulation établit des dépendances pour l'évolution de chaque espèce, spécifiées par les arêtes dirigées du graphe. Les dépendances seules, cependant, ne suffisent pas pour spécifier la fonction logique régissant l'évolution d'une espèce. Au lieu de cela, les fonctions logiques associées à chaque espèce, contraintes par le graphe d'influence, sont codées dans les paramètres d'un réseau de régulation discret. L'espace des fonctions logiques admissibles est alors représenté par un réseau de régulation paramétrique. D'une part, les réseaux de régulation paramétriques peuvent être utilisés pour l'identification de valeurs de paramètres pour lesquelles le réseau de régulation discret résultant satisfait des propriétés (dynamiques) données. L'identification des paramètres des réseaux de régulation peut ainsi être vue comme un exemple particulier de synthèse de modèle, dans le cadre contraint du graphe d'influence sous-jacent. D'autre part, les réseaux de régulation paramétriques peuvent être analysés comme un modèle autonome, pour faire des prédictions robustes vis-à-vis de la variabilité du réseau.L'analyse de la dynamique du réseau de régulation paramétrique est entravée par la double explosion combinatoire, de l'espace d'états et de l'espace des paramètres. Dans cette thèse, nous développons de nouvelles méthodes d'analyse de réseau de régulation paramétrique, sous forme de sémantique spécialisée, visant à atténuer l'explosion combinatoire. Tout d'abord, nous introduisons une interprétation abstraite de l'ensemble des évaluations de paramètres admissibles (paramétrisations). L'abstraction permet de représenter n'importe quel ensemble de paramétrisations par un encodage de taille constante, au prix d'une sur-approximation conservatrice. Deuxièmement, nous élevons la sémantique d'ordre partiel sous la forme d'un déploiement des réseaux de Petri vers des réseaux de régulation paramétriques. Les graphiques d'influence des réseaux de régulation biologique ont tendance à être relativement clairsemés, ce qui permet une grande concurrence. Cela peut être exploité par des méthodes de réduction d'ordre partiel pour produire des représentations d'espace d'état concises.Les deux approches visent à aborder les deux aspects de la double explosion combinatoire et sont introduites de manière compatible, ce qui permet de les utiliser simultanément. Une telle application est soutenue par une implémentation prototype utilisée pour mener des expériences sur divers réseaux de régulation paramétriques. Nous considérons en outre des raffinements des méthodes, comme une méthode de réduction de modèle à la volée portée aux réseaux de régulation paramétriques à partir de réseaux d'automates

Most Permissive Semantics of Boolean Networks

As shown in (http://dx.doi.org/10.1101/2020.03.22.998377), the usual update modes of Boolean networks (BNs), including synchronous and (generalized) asynchronous, fail to capture behaviors introduced by multivalued refinements. Thus, update modes do not allow a correct abstract reasoning on dynamics of biological systems, as they may lead to reject valid BN models.This technical report lists the main definitions and properties of the most permissive semantics of BNs introduced in http://dx.doi.org/10.1101/2020.03.22.998377. This semantics meets with a correct abstraction of any multivalued refinements, with any update mode. It subsumes all the usual updating modes, while enabling new behaviors achievable by more concrete models. Moreover, it appears that classical dynamical analyzes of reachability and attractors have a simpler computational complexity:- reachability can be assessed in a polynomial number of iterations. The computation of iterations is in NP in the very general case, and is linear when local functions are monotonic, or with some usual representations of functions of BNs (binary decision diagrams, Petri nets, automata networks, etc.). Thus, reachability is in P with locally-monotonic BNs, and P$^{\text{NP}}$ otherwise (instead of being PSPACE-complete with update modes);- deciding wherever a configuration belongs to an attractor is in coNP with locally-monotonic BNs, and coNP$^{\text{coNP}}$ otherwise (instead of PSPACE-complete with update modes).Furthermore, we demonstrate that the semantics completely captures any behavior achievable with any multilevel or ODE refinement of the BN; and the semantics is minimal with respect to this model refinement criteria: to any most permissive trajectory, there exists a multilevel refinement of the BN which can reproduce it.In brief, the most permissive semantics of BNs enables a correct abstract reasoning on dynamics of BNs, with a greater tractability than previously introduced update modes

Most Permissive Semantics of Boolean Networks

As shown in (http://dx.doi.org/10.1101/2020.03.22.998377), the usual update modes of Boolean networks (BNs), including synchronous and (generalized) asynchronous, fail to capture behaviors introduced by multivalued refinements. Thus, update modes do not allow a correct abstract reasoning on dynamics of biological systems, as they may lead to reject valid BN models.This technical report lists the main definitions and properties of the most permissive semantics of BNs introduced in http://dx.doi.org/10.1101/2020.03.22.998377. This semantics meets with a correct abstraction of any multivalued refinements, with any update mode. It subsumes all the usual updating modes, while enabling new behaviors achievable by more concrete models. Moreover, it appears that classical dynamical analyzes of reachability and attractors have a simpler computational complexity:- reachability can be assessed in a polynomial number of iterations. The computation of iterations is in NP in the very general case, and is linear when local functions are monotonic, or with some usual representations of functions of BNs (binary decision diagrams, Petri nets, automata networks, etc.). Thus, reachability is in P with locally-monotonic BNs, and P$^{\text{NP}}$ otherwise (instead of being PSPACE-complete with update modes);- deciding wherever a configuration belongs to an attractor is in coNP with locally-monotonic BNs, and coNP$^{\text{coNP}}$ otherwise (instead of PSPACE-complete with update modes).Furthermore, we demonstrate that the semantics completely captures any behavior achievable with any multilevel or ODE refinement of the BN; and the semantics is minimal with respect to this model refinement criteria: to any most permissive trajectory, there exists a multilevel refinement of the BN which can reproduce it.In brief, the most permissive semantics of BNs enables a correct abstract reasoning on dynamics of BNs, with a greater tractability than previously introduced update modes

Parameter Space Abstraction and Unfolding Semantics of Discrete Regulatory Networks

International audienceThe modelling of discrete regulatory networks combines a graph specifying the pairwise influences between the variables of the system, and a parametrisation from which can be derived a discrete transition system. Given the influence graph only, the exploration of admissible parametrisations and the behaviours they enable is computationally demanding due to the combinatorial explosions of both parametrisation and reachable state space. This article introduces an abstraction of the parametrisation space and its refinement to account for the existence of given transitions, and for constraints on the sign and observability of influences. The abstraction uses a convex sub-lattice containing the concrete parametrisation space specified by its infimum and supremum parametrisations. It is shown that the computed abstractions are optimal, i.e., no smaller convex sublattice exists. Although the abstraction may introduce over-approximation, it has been proven to be conservative with respect to reachability of states. Then, an unfolding semantics for Parametric Regulatory Networks is defined, taking advantage of concurrency between transitions to provide a compact representation of reachable transitions. A prototype implementation is provided: it has been applied to several examples of Boolean and multi-valued networks, showing its tractability for networks with numerous components

Unfolding of Parametric Boolean Networks

International audienceIn systems biology, models of cellular regulatory processes such as gene regulatory networks or signalling pathways are crucial to understanding the behaviour of living cells. Available biological data are however often insufficient for full model specification. In this paper, we focus on partially specified models where the missing information is abstracted in the form of parameters. We introduce a novel approach to analysis of parametric logical regulatory networks addressing both sources of combinatoric explosion native to the model. First, we introduce a new compact representation of admissible parameters using Boolean lattices. Then, we define the unfolding of parametric Boolean networks. The resulting structure provides a partial-order reduction of concurrent transitions, and factorises the common transitions among the concrete models. A comparison is performed against state-of-the-art approaches to parametric model analysis